CHANDRA observation of the narrow-line Seyfert 1 galaxy 
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Abstract 

The narrow-line Seyfert 1 galaxy (NLS1) IRAS 13224-3809 has been observed with the CHANDRA 
High Resolution Camera (HRC-I) for 12 ksec on February, 2, 2000. The source was proposed for CHAN- 
DRA observations to precisely determine the X-ray centroid position and to investigate the timing prop- 
erties of the most X-ray variable Seyfert galaxy. The position derived from CHANDRA confirms that the 
X-ray emission is associated with IRAS 13224-3809. The CHANDRA HRC-I light curve shows indications 
for a possible presence of a quasi-periodic oscillation. The strongest signal is found at 2500 sec. Accretion 
disk instabilities may provide a plausible explanation for the quasi-periodic oscillations. 

KEYWORDS: galaxies: active — galaxies: individual: IRAS 13224-3809 — galaxies: Seyfert — X-rays: 
galaxies 



1. Introduction 

The narrow-line Seyfert 1 galaxy (NLS1) IRAS 13224- 
3809 was the first galaxy proposed for monitoring obser- 
vations with the ROSAT High Resolution Imager (HRI) 
(Boiler et al., 1997). The most important observational 
fact is the detection of multiple strong flaring events dur- 
ing the 30 day observations. The authors found five 
giant-amplitude variations with the most extreme vari- 
ability of a factor of 57 in 2 days. In addition, other 
ROSAT observations were performed, mostly with the 
HRI detector during 1992 and 1998. To measure the 
precise X-ray position, IRAS 13224-3809 has been pro- 
posed for an observation with the CHANDRA satellite. 



2. Data reduction 

IRAS 13224-3809 has been observed with the Chandra 
High Resolution Camera (HRC-I) on February, 2, 2000 
between 15:16:28.59 UT and 18:41:04.25 UT with an ef- 
fective exposure time of 12271.6 sec. 

Inspection of the image with the source con- 
structed from the reprocessed level 2 photon event file 
(hrcf00328N002_evt2.fits) indicates that the aspect re- 
construction was not yet perfect as the FWHM of the 
source is ~ 1.04 arcsec which is larger than that expected 
from the instrument. A detailed analysis revealed that 
the source moved around the mean source position by 
approximately ± 0.5 arcsec with the dither period. This 
analysis was performed by computing the centroid in x, 



y, t for 100 consecutive events that lie within a circle 
of 2 arcsec around the mean source position. This 100 
event window is shifted always by one event giving a 
smooth curve which shows the displacement of the cen- 
troid versus time giving us a correction vector. This 
vector is then used to correct the position of the indi- 
vidual photon events. The data reduction, correction, 
and extraction was performed using routines written in 
IDL. Using the correction vector we obtain a more re- 
alistic FWHM of - 0.64 arcsec for IRAS 13224-3809. 
This value is still larger than the expected value for the 
instrument indicating that extended X-ray emission is 
present around IRAS 13224-3809. A maximum likely- 
hood source detection usinig CELLDETECT (from the 
CIAO data analysis package) performed on the binned 
image also indicates that the source FWHM is greater 
than that of the instrumental PSF (PSF-ratio=1.43187). 

After checking the radial profile of the source we chose 
an extraction radius of 9 arcsec around the source so 
that also the photons from the diffuse emission could 
be collected. The background selection was performed 
by selecting photons in 8 regions, each with a radius of 
9 arcsec. The regions are located 20 arcsec from the 
source. 

During the observation the total number of events, es- 
pecially all background, rose by a factor ~ 3-4 over a 
time span of ~ 2000 sec. Higher than normal solar ac- 
tivity is the most probable cause for this increase, which 
leads to a leakage of source counts based on the dead- 
time of the instrument. The ratio of detected to lost 



photons versus time is given in deadtime correction file. 
This file was then used to correct lower count rate of the 
source during times with high background count rates. 
Only the photons lying in the goodtime intervals were 
used. 

3. The CHANDRA position of IRAS 13224-3809 

With the CHANDRA satellite one can determine the 
position of the X-ray emission with an unprecedented 
precision. In Fig. [l], we show an optical image taken 
from the Palomar Digitized Sky Survey (DSS) overlaid 
with a CHANDRA X-ray contour map of IRAS 13224- 
3809. 
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Right Asesnsion (J2000) 

Fig. 1. Optical image (taken from DSS) overlaid with CHANDRA 
X-ray contour map of IRAS 13224-3809 and marked CHAN- 
DRA X-ray position of the sources GSC 7787 01660 (II) and 
RX J132537-3825 (III). 

The CHANDRA X-ray positions of two other sources, 
GSC 7787 01660 (II) and RX J132537-3825 (III), are 
marked by a cross. Note that no optical counterpart is 
found in optical catalogues for the source RX J132537- 
3825 (III). For the source detection we have used 
CELLDETECT from the CHANDRA CIAO data analy- 
sis package. The centroid position for IRAS 13224-3809 
in the CHANDRA image, computed with CELLDE- 
TECT, is a 20 oo = 13 fl 25T9.41", <5 200 o = -38°24'52.85". 
The coordinates of IRAS 13224-3809 and the two other 
X-ray sources are listed in table |l|. Fig. || shows an 
higher resolution optical image with the CHANDRA pos- 
tion marked by a red cross. The agreement between the 
X-ray positions and optical positions obtained from the 
GSC catalogue is within 2 arcsec. 




Fig. 2. CHANDRA X-ray positon marked with red cross in a deeper 
optical image (subrange of Fig. hi) with higher resolution observed 
at the ESO-Danish 1.54m telescope. 

4. Timing properties of IRAS 13224-3809 

In this section we discuss the properties of the light 
curve of IRAS 13224-3809, including the indications for 
a quasi-periodic signal and the tests for periodicity. 
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Fig. 3. Light curve of IRAS 13224-3809 combined with deadtime 
correction file (dtf) and goodtime correction file of the second 
processed data (correction files are shifted down by factor of i). 
Note the larger error in the_2000 sec region of dtf correction - 
for description see section bl The light curve indicate a quasi 
periodicity with a period of 2500 s. 

In Fig. I we show the CHANDRA light curve of 
IRAS 13224-3809. The reason for the lack of flaring 
events compared to the ROSAT observations in 1996 
(Boiler et al. 1997) is probably due the short expo- 
sure time (12 ksec) of the observation. However, the 
light curve exhibits flux variations with the strongest 
flux change by a factor 2 within 600 sec. The effi- 
ciency calculation shows that the change of luminosity 



Table 1. Coordinates: X-ray (using CELLDETECT) and optical coordinates (taken from HST GSC catalogue) of 
IRAS 13224-3809, GSC 7787 01660 (RX J132459-3826) and RX J132537-3825. 
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over time (AL/At = 3.2 • 10 41 ergs -2 ) exceeds the max- 
imum allowed value by accretion onto a Schwarzschild 
black hole. This fact points to the presence of a Kerr 
black hole in IRAS 13224-3809. The 'peak emission' of 
the light curve is not well defined and flattened, however 
four well-defined count rate minima are clearly visible 
and separated by 2500 sec. This fact and the shape of 
the count rate variations seems to be an indication for 
a quasi-periodic signal in IRAS 13224-3809. In Fig. |, 
the light curve with an inverted count rate axis illustrates 
the sharp and clearly separated peaks at 900, 3400, 6000, 
8700 and 10800 seconds. 



F x = (4.310 ±0.035) 
to the X-ray luminosity Lx 



10 n crg cm 2 s 1 



corresponding 
(3.880±0.036)-10 44 erg s" 1 



(for z = 0.0667, q = \ and H = 75 km s^Mpc -1 ). 

4.2. Tests of periodicity - x 2 , FFT, Lomb - Scargle 

To probe the quasi-periodic oscillation we used four in- 
dependent tests. In Fig. |^, the strongest signal in 
the Lomb-Scargle periodogram is 2500 sec. The second 
strong signal at 1300 sec seems to be the higher order of 
the signal at 2500 sec. The same result is given by the 
power spectrum of the source in Fig. o. 
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Fig. 4. Light curve with inverse count rate axis of IRAS 13224—3809 
illustrate the sharp and clearly separated peaks. The light curve 
indicate a quasi periodicity with a period of 2500 s. 

The mean count rate of IRAS 13224-3809, obtained 
from the light curve, is cps = 1.222 ± 0.010 counts s _1 . 
The unabsorbed flux in the 0.1 - 10.0 keV energy range 
computed with W3PIMMSQ using a power law and black 
body model with the spectral parameters; hydrogen col- 
umn density N H = (0.87 ± 0.05) • 10 21 cm~ 2 (Boiler et 
al., 1996), photon index T = 2.47, temperature kT = 
0.121 keV (Brandt, private communication) and the rel- 
ative normalisation n = 0.2181 at 1.0 keV; results in 

* 1 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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Fig. 5. Lomb scargle test for periodicity in the 0.2-10 keV energy 
range. 

We have folded the light curve with different periods, 
ranging from 500 to 3500 seconds, and have determined 
the corresponding x 2 value. The reduced x 2 versus the 
folding period in Fig. [7| shows also the strongest peak at 
2500 sec. The corresponding folded light curve is given 
in Fig. U with the best fitting period of 2400 sec and 
modulations of about 20 %. We made simulations of 
light curves with the observed time sequence and phase 
randomized for a red noise f~ l power spectrum (see Fig. 
|9|). The indication for the periodic signal in the light 
curve of IRAS 13224-3809 (see the power spectrum in 
Fig. o) arc in the same order compared to the power 



seen in the simulated power spectrum. Therefore we 
need the 80 ksec XMM-Newton observation of the source 
to confirm the presence of any quasi-peridic signal. 
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Fig. 8. The light curve of IRAS 13224-3809, folded with a period of 
2380 s. 



Fig. 6. Power spectrum for IRAS 13224-3809 in the 0.2-10 keV 
energy range. 
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Fig. 7. Reduced \ 2 versus the folding period for the 0.2-10 keV 
energy range. 



5. Discussion 

The quasi-periodic oscillations if confirmed by other X- 
ray observations, can be caused by instabilities of the 
inner parts of the accretion disk, which can lead to vari- 
ations of the accretion rate. In this case the period is 
expected to be of the order of the radial drift time scale 
(" instability time scale" ) 

a \RsJ c 

with e ~ 5 — 10 and a ~ 0.1 (Sunyaev, private com- 
munication). For a period of 2500 sec and a mass of 
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Fig. 9. Simulated red noise power spectrum of IRAS 13224—3809, 
normalized to the real data set. 



M BH d ~ 3 • 10 6 Mq derived from the Eddington limit, 
we have estimated a radius of r w 1.4 Rg, where the in- 
stabilities occur, which would be also indicative for the 
presence of a Kerr black hole. Absorption by Compton- 
thick matter can probably ruled out by an estimation for 
the distance of the absorbing matter. 

In Fig. [l^ we show the lower and upper value of the 
black hole mass for IRAS 13224-3809. The lower mass 
is given by the Eddington limit (M^ cr ~ 3 ■ 10 6 M Q ). 
The upper mass limit is estimated from material orbit- 
ing with v = c ( MBPj por - 4 • 10 7 M ). The blue curve 
demonstrates the dependence of black hole mass to dif- 
ferent rotation radii in units of the Schwarzschild radius, 
assuming that rotating Compton-thick matter is causing 
the X-ray periodicity (semi-relativistic calculation)^. 
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r in units of Rg and v = 



27TI . = . 
T ' 



Rs 



The putative Compton-thick absorbing material has 
to be located at distances of ~ 5 Rs- This seems to 
be inconsistent with neutral Compton-thick absorbing 
material, because it is too close to the central black hole. 



Boiler Th., Brandt W.N., Fabian A.C., Fink H.H., 1997 

MNRAS, 289, 393 
Sunyaev R.A., 1973, Soviet Astron. AJ, 16, 941 
Sunyaev RA., 2000, private communication 




Fig. 10. Dependence of black hole mass to different rotation radii for 
rotating Compton-thick matter (blue curve). The radius is given 
in units of the Schwarzschild radius and the black hole mass in 
units of solar mass. 



6. Summary 

The precisely determined X-ray positon of the flaring 
source IRAS 13224-3809 (Boiler et al. 1997) is in good 
agreement with the optical position from the Hubble 
GSC of the galaxy (within 2 arcsec). The light curve 
of IRAS 13224-3809 shows no huge flaring events as 
seen by Boiler et al. (1997); however, we found indi- 
cations for a quasi-periodic signal in the light curve of 
IRAS 13224-3809 with a period of 2500 sec. This quasi- 
periodic signal might be caused by instabilities of inner 
parts of the accretion disk. To confirm and to study 
the triggering processes of the putative periodic signal 
we need the scheduled XMM-Newton observation which 
will yield both timing and spectral information of this 
source. 

In addition, we found evidence for the presence of a 
Kerr black hole in IRAS 13224-3809. The first hint is 
given by the change of luminosity over time, which ex- 
ceeds the maximum allowed value by accretion onto a 
Schwarzschild black hole. The second one is the estima- 
tion of the distance of the accretion disk instabilities to 
the black hole, resulting in the short distance of about 
1.4 Schwarzschild radii. We estimate the black hole mass 
of IRAS 13224-3809 between 3 • 10 6 M„ and 4 • 10 7 M fl . 
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